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In generic particle physics models, the inflaton field is coupled to other bosonic and fermionic fields 
that acquire large masses during inflation and may decay into light degrees of freedom. This leads 
to dissipative effects that modify the inflationary dynamics and generate a nearly-thermal radiation 
bath, such that inflation occurs in a warm rather than supercooled environment. In this work, we 
perform a numerical computation and obtain expressions for the associated dissipation coefficient, 
focusing on the regime where the radiation temperature is below the heavy mass threshold, while 
allowing for generic couplings and field multiplicities. We also discuss radiative corrections to particle 
masses and decay widths and how they may significantly increase the effects of dissipation. 



I. INTRODUCTION 



Inflation [T| is an extremely successful paradigm providing an elegant solution to the shortcomings of the standard 
cosmological model, in particular explaining the high degree of flatness and homogeneity of the observable universe, 
as well as describing the origin of the temperature anisotropies observed in the Cosmic Microwave Background and 
the seeds of the Large Scale Structure of our universe. 

In the traditional picture of inflation, the early universe is dominated by the vacuum energy of a scalar field <f> 
that slowly rolls down its potential, V ((/)). This mimics the effect of a cosmological constant if the kinetic energy of 
the inflaton field is negligible, a condition that is necessarily violated at some point during the evolution of the field 
towards the minimum of the potential, thus yielding a finite period of accelerated expansion. This typically requires 
a scalar field with extremely weak self-interactions, so that its very flat potential leads to 40-60 e-folds of inflation. 
The quasi-exponential expansion redshifts away any matter or radiation present before the scalar field comes to 
*^ 1 dominate the energy density: 

5 : p R + 4Hp R = => PR ~e- iHt , (1.1) 

and consequently one needs a mechanism to reheat the cold inflationary universe and provide an exit into the standard 



radiation-dominated cosmology. The inflaton field must thus have some interactions with other fields, which directly or 
indirectly lead to the production of at least the Standard Model degrees of freedom and the particles that presumably 
constitute the inferred dark matter content of our universe. 

These interactions must, however, be sufficiently weak in order to preserve the flatness of the inflaton potential at 
k> ■ the quantum level, which can be achieved by either (fine-)tuning the associated couplings or by considering additional 
; 1 1 symmetries that, in particular, keep the mass of the scalar inflaton below the value of the Hubble parameter during 
C$ , inflation, thus ensuring an overdamped evolution. Moreover, when one considers generic couplings to additional scalar 
(fermionic) fields of the form g 2 4> 2 x 2 (s^VOj these fields typically become heavy during inflation due to the large 
values of the inflaton vacuum expectation value (vev). This implies that the decays <fr —> xXjV'V' are generically 
forbidden in the slow-roll regime, with reheating of these degrees of freedom occuring naturally only at the end of 
inflation if the inflaton vev drops sufficiently. 

Interactions with additional scalar fields may also have interesting dynamical consequences, as is the case of hybrid 
inflation models @, H|, where negative contributions to the squared mass of the waterfall field x, coupled to the 
inflaton as above, make it tachyonic at some critical value of <j). Inflation then ends with a phase transition in which 
the x -ne ld evolves towards the true minimum of its potential, and in the process it may reheat the universe if it is 
coupled to radiation. 
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Warm inflation [J, H| (see also @) provides an alternative picture of the inflationary universe, in which particle 
production is sourced by the rolling inflaton field itself, so that radiation is not completely diluted away: 

T6 2 

p R + AH p R = T4> 2 => PR^-rrr- const - > (1-2) 

in the slow-roll regime. Although accelerated expansion only occurs for a sub-dominant radiation component, pn <C 
Pef,, it is possible that its temperature exceeds the de Sitter temperature, T > H, which significantly changes the 
inflationary dynamics 1 . On one hand, the spectrum of primordial density fluctuations is seeded by thermal rather 
than vacuum fluctuations of the inflaton field [J, I7l4l0| , which my lead to interesting observational consequences such 
as the suppression of the tensor-to-scalar ratio and significant deviations from a gaussian spectrum [Ilr-|l4f . On the 
other hand, particle production induces an additional friction term in the inflaton's motion: 

P4, + 3H(j>t+ P4 ,) = -T4> 2 , (1.3) 

which using p^ = <j) 2 /2 + V{<f>) and = <j) 2 /2 — V(<f>) leads to: 

4> + 3H^ + V ct> = -T<j> , (1.4) 

where denotes the derivative of the inflaton potential. This friction then contributes to overdamp the inflaton's 
motion, alleviating the need for very finely-tuned flat potentials (see e.g. [3]). This may be particularly relevant 
for supergravity and string theory models fl6l where one typically finds > H , precluding a sufficiently long 
period of inflation. 

We are mainly interested in dissipative effects in the adiabatic regime, where the microscopic particle dynamics 
is much faster than the evolution of any macroscopic variable, which typically holds during slow-roll inflation. The 
calculation of dissipation coefficients in the adiabatic regime has been a long studied problem in quantum field theory, 
starting primarily with the seminal works in the 80's of Hosoya and Sakagami [lf| for the cj) 4 interaction (see also [l9j). 
This was followed by Morikawa [2(|, who used the Closed Time Path formalism and obtained an effective Langevin- 
like equation, including an explicit fluctuation-dissipation relation. Fluctuation-dissi pat ion relations emerging from 
quantum field theory models have since been examined by several other authors |2l| - t24l |. Initially these treatments 
considered weak dissipation, leading to an underdamped evolution, although it was subsequently shown in (2f| that 
strong dissipation with overdamped trajectories could also be achieved. All these studies were performed in Minkowski 
spacetime and one of the first to consider dissipation in curved spacetime was Ringwald p6| , with other subsequent 
treatments such as (27l - [29j . Dissipation from quantum field theory models has in fact been considered for a variety 
of applications beyond warm inflation dynamics, such as phase transitions, heavy ion collisions and conventional 
reheating after inflation. 

For the reasons described above, the interactions between the inflaton and the radiation must be sufficiently sup- 
pressed to maintain the flatness of the potential, and in particular the finite temperature of the radiation may induce 
large thermal corrections to the inflaton mass, with m$ ~ T > H . This in fact precludes successful realizations of 
warm inflation in the simplest models, with the inflaton coupled directly to light scalars or fermions as described 
above, as the additional friction cannot overcome the increase in the inflaton's mass [25|,[3(|. This of course assumes 
these are relativistic degrees of freedom, with T 3> rn x ^, but, as previously discussed, fields which are directly coupled 
to the inflaton tend to acquire large masses during inflation. 

This has motivated implementations of warm inflation scenarios in the low-temperature regime, with m x ^ < T (see 
e.g. [3l[ for other recent implementations of warm inflation). In this case the on-shell production of these degrees of 
freedom is Boltzmann-suppressed, but contributions from virtual modes may still lead to significant dissipative effects 
if the fields have a finite decay width. In fact, this scenario, proposed in [32| and suggestively known as the two-stage 
mechanism, has all the ingredients of hybrid inflation models, where it is the waterfall field and not the inflaton that 
interacts with and may decay into light degrees of freedom. Moreover, it has been shown that the inclusion of both 
bosonic and fermionic degrees of freedom in a supersymmetric theory may help controlling both radiative and thermal 
corrections to the inflaton potential, despite supersymmetry being broken by both the finite energy density and the 
finite temperature during warm inflation (33j . 

In this work, we extend earlier computations of the dissipation coefficient T in the low-temperature regime (34l . |35| 
for the case where the fields coupled to the inflaton have a large decay width, where one expects particle production 
to be more efficient. Although most computations have to be performed numerically in this regime, our goal is to 



1 In this work we will consider a radiation bath that is close to thermal equilibrium, although this need not be the case in general. 
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provide accurate expressions for T that may be used in constructing models of warm inflation for a broad range of 
field masses and couplings. 

Our work is organized as follows. In the next section we describe the interactions between the inflaton field, the 
heavy fields and the light degrees of freedom that make up the radiation bath. We then discuss the leading thermal 
corrections to the particle masses for weak and strong couplings, which may be found in Section III. In Section IV, we 
describe the computation of the dissipation coefficient in the low-temperature regime and determine its dependence 
on the particle masses and couplings by fitting our numerical results to well-motivated expressions. Finally, in Section 
V we estimate the effects of higher-order radiative corrections to the decay width of the heavy fields and discuss their 
implications for warm inflation model building. We summarize our main results and conclusions in Section VI. 



II. SUPERSYMMETRIC MODEL 



We consider a generic supersymmetric model with chiral superfields $, X and Y i7 i = 1, . . . ,Ny, described by the 
superpotential [HI, HH : 

W = U.X 2 + *±XY? + /($) , (2.1) 

where a sum over the index i is implicit. The scalar component of the superfield $ describes the inflaton field, with 
an expectation value <f> = ip/y/2, which we assume to be real, and the generic holomorphic function /($) describes the 
self-interactions in the inflaton sector. The structure of the superpotential is quite generic in rcnormalizablc models, 
following from the simple assumption that not all degrees of freedom are directly coupled to the inflaton field, and 
may for example be naturally implemented in D-brane constructions and related gauge theories (see e.g. [f7j and 
references therein). The analysis that we will go through below concentrates mostly on the scalar and fermionic 
sectors, which are already sufficiently involved, but we discuss the possibility of embedding this superpotential in a 
gauge theory in Section VI. 

The Lagrangian density describing the interactions between the inflaton vev and the scalar components of the 
superfields X and Yi, denoted by x an d respectively, is given by: 

Scalar = V(<p) + \ 9 V \X? + f ( X 2 + X *) + ^ | X | 4 + 

+ t5 H 2 + xtff2 ) + hJ r^ + h ^ 2 ^ 2 ' (2 - 2) 

where V{ip) — |/'(<?!>)| 2 is the potential driving inflation. Similarly, the interactions involving the fermionic components 
ip x and ip ai and the scalar inflaton are given by: 

Zfermion = ^ P L *1>X + yXVV; Pl VV; + ^CTiVv, Pl lp X + n ' C ' ' ( 2 - 3 ) 

where Pl = (I — 75) /2 is the left-handed chiral projector. Note that for our study of the dissipative dynamics during 
inflation we are only interested in the interactions involving the inflaton vev, but one should take into account that 
there are also interactions involving scalar fluctuations about the background value and their fermionic superpartners, 
although, for simplicity, we do not write them explicitly in Eqs. (|2.2I) and (|2.3[) . Also, for the purposes of this work, 
we will assume that the couplings g and hi are real, although complex couplings may play an important role in, for 
example, generating a baryon asymmetry through dissipative effects during warm inflation [36j . 

From the interactions shown in Eqs. (|2.2p and (|2.3p . we can see that the fields in the X multiplet, which are 
directly coupled to the inflaton in the superpotential, acquire tree-level masses that are generically large for large 
values of the inflaton vev. On the other hand, the fields in the Yi multiplets remain massless at tree-level, although we 
may in general consider an additional mass term for these fields, which we will assume to be small compared to the 
temperature scale during inflation, an issue we will return to in the next section. The non- vanishing potential energy 
driving inflation breaks supcrsymmetry, which leads to a mass splitting in the scalar sector of the X multiplet: 



lxR '< 2 y gtp 2 J ' 

= ^— , (2.4) 
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where x ~ (xr + *X/)/v^2. Due to the smallness of the inflaton self-interactions typically required for a sufficiently 
long period of slow-roll inflation and a small enough amplitude for the spectrum of primordial density perturbations, 
this splitting may in general be neglected during inflation. For example, for V(ip) = \ 2 (p 4 , the splitting factor is 
1 ± X/g, with the primordial spectrum imposing A ~ 10~ 7 , so that we generically expect A <C g (see e.g. [HI). 
Henceforth we will then consider a common mass m x = gip/^/2 for all fields in the X multiplet. If the inflaton vev is 
sufficiently large, this mass will be above the temperature of the radiation bath during warm inflation, in which case 
thermal contributions to the inflaton mass from \ an d ip x loops are Boltzmann-suppressed, e~ m */ T <C 1, and do not 
destroy the flatness of the tree-level potential. 

This mass splitting may nevertheless be important when computing radiative corrections to the inflaton potential, 
and the supcrsymmctric Coleman- Weinberg potential has been shown to yield 



Vcw = TTT~2 Str 



M\ log 



32ir 




' V(<p)\og[ ^ I , 



where fi is the renormalization scale and we have taken the leading correction from supersymmetry breaking. This 
shows that the leading radiative corrections to the inflaton potential are logarithmic and hence will not spoil the 
flatness of the potential for g < 1. Note that if we consider Nx multiplcts coupled to the inflaton as in Eq. (12.11) . 
the radiative corrections will be proportional to g 2 Nx, and this may still accommodate a moderately large field 
multiplicity without a significant tuning of the coupling constant. 

One should also notice that in some cases the mass splitting induced by supersymmetry breaking may actually 
become relevant at the end of the inflationary evolution. For example, for /(<£>) = jM 2 $, we have V(ip) ~ g 2 M 4 and 
m XR i = m x -^ 2 /2, so that the imaginary component becomes tachyonic for ip < AI, leading to a phase transition 
that ends inflation in the conventional hybrid mechanism. 

The superpotential (|2.1j) is thus a simple generalization of supcrsymmetric hybrid inflation that allows for an 
arbitrary inflaton potential and additional couplings of the waterfall ficld(s) to other scalars and fermions. These 
couplings allow for the decays x <J<J ^ ^Pa^Pa and ip x —> cripo-, as can be easily seen from Eqs. (|2.2[) and (|2.3[) . and thus 
allow the waterfall field and its superpartners to 'reheat' the universe after the hybrid transition. Most importantly, 
even though the fields in the X-multiplet are too heavy to be produced on-shell during inflation, virtual x an d ipx 
pairs can be produced from the rolling inflaton field and may in turn decay into the F^-multiplef particles, leading 
to dissipative particle production and to the two-stage realization of warm inflation [32j , which will be main topic of 
this work. 

The scalar and fcrmionic Lagrangian densities also include interactions involving the fields in the X and Yi multiplcts 
that play an important role in the dynamics of warm inflation, not only by contributing to the thermal masses of 
these fields, as we will analyze in the next section, but also by allowing thermal scatterings that keep the system close 
to thermal equilibrium. For example, the cross-section for thermal scatterings of the complex light scalars <ii includes 
contributions from UiUi — > <Jj(Jj and Ui<j\ — > cr^ crj , being given by: 

< - *3£±2 , (2-6) 

where y/s is the centre-of-mass energy and we have assumed hi = h for all species, which for simplicity we will consider 
in the remainder of our discussion. For ^Js ~ T during warm inflation, this yields a scattering rate: 

3..te4^ m+J) (|) Km (*)V(|), 

for Ny ^> 1, so that we expect T\ h > H for T > H and perturbative couplings. On the other hand, fermionic 
species can only interact with each other or their scalar superpartners via the exchange of virtual x, ^ x -modes, with 
scattering rates thus suppressed by a factor (T/m x ) 4 <C 1. As we shall discuss below, fcrmionic contributions to the 
dissipation coefficient are also suppressed by similar factors, so that we expect fermions to be a negligible component 
of the radiation bath, a symptom of supersymmetry breaking during warm inflation. 



III. THERMAL MASS CORRECTIONS 



The two-stage interactions between the inflaton and X and Yi fields leads to the formation of a thermal bath 
during inflation, and finite temperature effects may thus become important. Since x and ip x acquire large tree- level 
masses and we assume T m x , their contributions to thermal loop corrections are Boltzmann-suppressed and may 
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be ignored. This precludes large thermal mass corrections not only to the inflaton field but also to the light a and ip a 
particles that form the radiation bath. 

The masses of both the bosonic and fcrmionic components of the X-multiplct receive nevertheless corrections from 
their interactions with the light fields, the leading diagrams being illustrated in Fig. [TJ 




FIG. 1: Feynman diagrams contributing to the self-energy of the ip x and \ fields to leading order. Note that the last diagram 
corresponds to an effective finite temperature contribution to the coupling g(T) rather than a thermal mass correction [331 ] . 

These have been computed in [33j and are identical for both bosonic and fermionic degrees of freedom. Taking into 
account our normalization of the couplings in the superpotential (which differs from [331 ] . by a factor 1/2) and the 
field multiplicity, this yields the renormalized mass: 

m 2 = m 2 + ^T 2 . (3.1) 

In particular, this implies that the one-loop effective potential at finite temperature can be obtained by replacing 
the treedevel mass by m x in Eq. (|2.5p . This shows that the leading thermal corrections to the inflaton potential are 
logarithmic and that, for T <C m x , they are smaller than the zero-temperature radiative corrections induced by the 
inflaton vev. We will nevertheless include these corrections in the computation of the dissipation coefficient. 

In the Yrsector, the light er.j scalars receive thermal corrections from their self-interactions, while their supcrpartners 
only have couplings involving at least one heavy field, and thermal loop contributions are thus suppressed in this case. 
The leading diagrams contributing to the scalar masses are shown in Fig. [2j 




FIG. 2: Feynman diagrams contributing to the self-energy of the Oi fields to leading order. Note that the 'sunset' diagram on 
the right gives the leading-order correction involving different light species. 

The leading corrections to the mass of each species cr, correspond to the quartic self-interactions (/i 2 /4)|<7j| 4 , 
while the mixing between different species, (hihj /4)fT 2 crj 2 gives rise to the higher-order 'sunset' diagram in Fig. [2j 
proportional to h A Ny 2 ■ To obtain the physical o^-mass one needs to solve the associated gap equation, and to lowest 



2 Note that the relevant two-point correlation function for a complex scalar is (<7jcrt), which receives distinct contributions from the 
quartic and bi-quadratic interaction terms. 
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order we have [13, IH| : 

/TOct\ 2 _ /r«o \ 2 rn 2 h 2 
\T J ~\T J I6ir 2 

where we have omitted the field index for simplicity, mo denotes its bare mass and fi is the renormalization scale, 
with the physical mass evaluated at /i = m a . For small couplings, the hard thermal loop approximation yields 
m 2 ~ m,Q + h 2 T 2 / '12. We have solved the gap equation numerically for different values of m /T, the results being 
illustrated in Fig. [3] 



log 



h 2 
2^2 



dx 



Jx 2 + (m a lT\ 2 P V x2 +( m «/TV 



(3.2) 




As can be seen in Fig. [3l the hard thermal loop approximation holds for h < 0(1), and for large couplings the mass 
tends to a constant value, which is well described by the expression: 

m a 7r 

~¥ ~ 2 - ln(m 2 /T 2 ) - (3 ' 3) 

This shows that, for mo % T, we expect m CT < T for all values of the coupling constant. This was expected due 
to the Boltzmann-factor in the loop integral, which precludes a large mass being generated by large couplings. Of 
course the one-loop approximation is not an accurate estimate for large couplings, but this result should hold at all 
orders in perturbation theory, including the 'sunset' contributions that may involve a large field multiplicity, as all 
loop diagrams involve this Boltzmann-supprcssion. Moreover, higher-order loop contributions do not have a definite 
sign (37j . which implies a slower growth of m a /T with the coupling constant and field multiplicities. 

The tree-level masses of the fields in the i^-multiplets have, of course, a model-dependent value, and may be 
generated by a variety of both supersymmctric and non-supcrsymmctric effects which are not directly related to the 
inflationary dynamics. However, if for example this sector includes the MSSM fields, we expect these masses to vanish 
or to at least be extremely small compared to the high temperatures typical of warm inflation models, T 3> O(TeV) 
[lj|. Given the above considerations, this will also be the case when thermal corrections are included, so that we will 
take m a < T in computing the associated dissipation coefficient. 



IV. DISSIPATION COEFFICIENT 



The coupling between the inflaton and the x, 4>x fields induces time non-local corrections to the inflaton effective 
action which, in the adiabatic regime, where fi/ip <SC T x , with the latter denoting the x decay width, lead to the 
effective friction term Tip in the inflaton's equation motion |35| . As mentioned above, these dissipative effects are due 
to the finite decay width of the X-multiplet fields, and correspond to an effective production of relativistic degrees of 
freedom, yielding a non-zero imaginary part for the self-energy of the fields in the l^-sector [39j . 
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The leading contribution to the dissipation coefficient from the complex scalar ^-modes arises at one-loop order, as 
illustrated in Fig. 21 and has the following form [351 ]: 



d 4 p 

W) 



A p^n B {l + n B ) 



(4.1) 



where ub{po) = [e p °/ T — 1] 1 is the Bose-Einstein distribution. 



/ \ 



g V i 
2 t 



FIG. 4: Feynman diagram contributing to the inflaton effective action at leading order. 



The spectral functions for the x-hclds correspond to the fully dressed propagators, including the effect of their finite 
decay width into Oi -particles: 



Px(Po,P) = 



4w p r x 



4uj 2 T 2 

p x 



(4.2) 



wher 



m\ 



p 2 for modes of 3- momentum |p| = p and energy po. The contribution from the i/^-modes is 

given by an analogous expression with fermionic propagators and the relevant couplings, but has been shown to be 
suppressed in the low-temperature regime, T <C m x [35j, which is the main scope of this work. It should be noted that 
the dissipation coefficient receives positive contributions from both the bosonic and fermionic fields that are directly 
coupled to the inflaton, so that only the time-local radiative corrections can be (partially) canceled by the underlying 
supersymmetry of the model. The suppression of the fermionic contribution is also a symptom of supersymmetry 
breaking during warm inflation, due to both the finite inflaton energy density and the finite temperature of the 
radiation bath. 

The decay width of the x bosons also includes contributions from both the bosonic and fermionic final states in 
the Fi-multiplets. The temperature and inflaton field dependence of the partial decay widths is also different for 
bosons and fermions (3o| . and again the fermionic contibutions are sub- leading. Taking also into account our earlier 
discussion of thermalization rates in the Y^-sector, it is clear that only the complex scalar fields \ an d o% are relevant 
to the dynamics of warm inflation, and that the underlying supersymmetry is only manifest in the logarithmic form 
of the radiative corrections to the inflaton potential in Eq. (|2.5[) . 

The leading process contributing to the decay width of the x fields is then the two-body decay x ~^ (and the 
complex conjugate process), and at finite temperature we include contributions from both decays and inverse decays, 
as well as thermal scatterings off particles in the radiation bath. This has been computed in 35[ from the imaginary 
part of the x _ self energy at one-loop order, yielding: 



h 2 N Y ml , 

64-7T LO n 



(4.3) 



where 



Ft(j>,Po) 



P 



T 
P 



log 



T 1 



1 -e 



T 1- er~ 1 - e 

-log — 

P \l-e-—l-e 



P0+" 



with 9{x) denoting the Heaviside function and 



PQ — " + 



{pl -p 2 - Ami) + 



{-pI+p 2 ) 



(4.4) 



P ± Po 1 



Ami 



Pl 



(4.5) 
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Note that the factor m 2 , appearing in Eq. (|4.3[) arises from the fact that the coupling between the x an d &i fields 
depends on the inflaton vev, whereas for the physical mass of the heavy bosons we use the renormalized mass rh x . 
Also, the first term in Eq. (|4.4j) corresponds to direct and inverse decays, while the second Landau-damping term is 
associated to thermal scatterings. 

From Eqs. (|4. Iti4.3|) . it is easy to see that the dissipation coefficient depends on the mass of \ fields, m x /T, and the 
effective coupling, h = hy/Ny, up to an overall rescaling by the coupling g 2 . We have integrated Eq. (|4.f j) numerically 
for different values of these parameters, and we show a representative selection of our results in Fig. [5j 
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FIG. 5: Numerical results for the dissipation coefficient as a function of the effective coupling h = hy/ Ny for different values 
of m x /T. These results correspond to m CT /T = 0.01. 

As one can easily conclude from the results shown in Fig. [5J the dependence on the effective coupling h is non-trivial, 
exhibiting a distinct behavior for small and large couplings depending on the mass of the \ fields. As we will discuss 
below, these correspond to two different contributions to the integral in Eq. (|4.f p that can be analyzed separately 
and, in fact, represent different physical processes. 



A. Low momentum contribution 



For m x /T > 10 and the range of values of the effective coupling shown in Fig. [5J the main contribution to the 
dissipation coefficient comes from virtual x modes with low momentum and energy, p,po <C m x , so that one can use 
the approximation (pq — lo 2 ) 2 ~ m x . If, in addition, these modes have a narrow width and thermal mass corrections 
can be neglected, T x <C rh x ~ m x , the spectral function takes the simple form p x ~ 4r x /m x and it is easy to see 
that: 

yLM / rp \ 3 

— =Ah i — . (4.6) 
9 m x \ m xJ 

The constant can be determined from the numerical data and one obtains A ~ 1.63 x 10 -3 for m a /T = 0.01, 
consistently with the results obtained in (34l . l3o| , taking into account our choice of normalization in Eq. (|2.I [I . This 
behavior corresponds to h < 1 for m x /T > 10. It is convenient to write this in the form: 

T LM = cJ-, , C = 0.00325/i 4 . (4.7) 

However, as one increases the effective coupling h, the finite width of the x fields and their thermally-induced mass 
become significant, and the dissipation coefficient grows slower than h A for h > 1 . Although the exact dependence 
cannot be determined analytically, the numerical data is well described by including a polynomial correction factor 
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in Eq. (J4T6J) , motivated by the r x -dcpendence of the loop integral in Eq. (|4.1[) and the form of the renormalized mass 



T 



LAI 



Ah 4 



9 z m x f(h,m x /T,m a /T) \m x 



T 



f(h,m x /T,m (T /T) = l + J2 



(4.8) 



where a n = ot n {m x /T, m a /T) and A = A{m rT /T). We have performed a numerical fit to the data using this expression 
in order to obtain the value and mass-dependence of these coefficients, by varying m x /T for fixed m a /T and vice- versa. 
The results for m a /T = 0.01 are illustrated in Fig. [SJ while the results for m x /T = 50 are shown in Fig. [7J 



io- 10 



m [r /T=0.01 











<*i ■ — . . , 











50 100 



500 1000 



FIG. 6: Coefficients of the polynomial correction to the low-momentum dissipation coefficient for different values of m x /T and 
rricr/T = 0.01. The solid lines represent the best fit curve in each case. 
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FIG. 7: Overall normalization and coefficients of the polynomial correction to the low-momentum dissipation coefficient for 
different values of m a jT and m x /T = 50. The solid lines represent the best fit curve in each case. 

From these results one can deduce the following approximate form for the polynomial coefficients and the overall 
normalization constant: 



-LM 



g 2 m x 



4.3 x 10" 4 ^ 4 | log(l - e-i^)| 



1 + 



OMh 2 + 0.02/t 4 (2?.) 2 + 0.002/16(1^) 4 + 7x 10-5/>(^)" 6 | |log(l-e- 2l ^)| 



— ) 



The expression (|4.9[) is valid for m x /T > 10, while for lower values one observes larger deviations from the best 
fit result, in particular for the lowest coefficient a 2l which exhibits a mild dependence on m x /T. This is due to 
contributions from real x-modes, which become significant for lighter x-hclds and smaller couplings and that we will 
analyze in the next section. 
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An important feature of these results is the mild logarithmic dependence on m a /T, and for practical purposes its 
is useful to evaluate this in the limit m a <C T, for which we obtain: 

T ^ * ^ Tl (4 10) 

1 + o.53^ (^f)~ 2 + 0.07^ 4 (^y A + 0.002^6 (^)~ 6 <^ 2 

In this limit and for m x /T > h, we recover an expression for the dissipation coefficient of the form in Eq. (|4.7[) : 

T LM ~ , C7 ~ 0.03/i 2 JV Y . (4.11) 



Eq. (|4.1ip is one of our main results, showing that, for large values of the effective coupling in the narrow width 
regime, the dissipation constant grows like oc h?Ny and not as h 4 N Y as obtained in [3_4j, l35( and which is valid 
only for h < 1. In Fig.[8]we also plot the effective function C<p(m x /T 1 m a /T) valid for m x /T > 10, illustrating how it 
becomes constant for m x /T > h and also showing that the dependence on m a /T is negligible for large values of the 
effective coupling. 





FIG. 8: The dissipation function C^(m x /T, m a /T) as a function of m x /T for different values of the effective coupling h — h\fNy 
and for m a /T — > (solid), m a /T = 0.01 (dashed) and ma/T — 1 (dotted), using the best fit expression in Eq. (|4.9|l . 



On the other hand, when m x /T < h, the effects of the width of the x fields become important and suppress the 
dissipation coefficient, as illustrated in Fig. [5] by the sharp decrease in C<f,. In the broad width limit, where the term 
with coefficient as is the leading polynomial correction, we have for m a <C T: 



-LM 



0.94 



9 6 ^ 

h 4 Ny T 3 



(4.12) 



B. Pole contribution 



The contribution to the dissipation coefficient from real x _m odes is generically Boltzmann-supprcssed, so that 
previous analysis assumed that it only becomes relevant in the high-temperature limit, T 3> m x . However, our 
numerical results suggest that the Boltzmann suppression factor e~ m ^l T may be compensated by a sufficiently small 
effective coupling in the low-temperature regime, so that on-shcll modes become the dominant contribution to the 
dissipation coefficient, as can be seen in Fig.[5]for m x /T < 10. To better understand this, let us expand the spectral 
function in Eq. (|4.2[) about its pole at p ~ lj p , which corresponds to the production of on-shell modes. We then 
obtain: 
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from which it is clear that T oc h~ 2 , a behavior that is evident in Fig. [5] in the small effective coupling region. 
Proceeding as in the previous section, we have performed a numerical fit to the data in this region to an expression 
of the form: 



— = B(m x /T,m cr /T)h' 

9 m x 



(4.14) 



The dependence on the mass of the light fields is actually negligible in this case for m a < T, and in Fig.|9]we illustrate 
our results for m a /T = 0.01. 

m o ./T = 0.01 



o.oi 



10 ' 



10 ' 



5 10 15 20 

FIG. 9: Value of B(m x /T) in the pole approximation for m a /T = 0.01. The solid red line corresponds to the fitting function 
B(x) = ie x /(e x - l) 2 . 

The coefficient B exhibits a strong dependence on m x /T, which is mainly due to the Bose-Einstein factors in 
Eq. (|4.ip . so that the dissipation coefficient can in this case be well approximated by: 



^£1 



-j 1 m x e mx ^ T 



I) 2 



(4.15) 



One should note that even though on-shell modes can give the dominant contribution to the dissipation coefficient 
for sufficiently small effective coupling, the width of the fields may be very small in this limit, which implies that the 
adiabatic condition (p/<p -C T x may no longer hold, which needs to be taken into account when constructing models of 
warm inflation in this limit. Also, it has been shown that a full resummation of ladder diagrams contributing at the 
same order to the dissipation coefficient is required in this limit, although this does not generically change its order 
of magnitude (4fjj |. 



C. Full dissipation coefficient 



While previously we have analyzed the contributions from virtual and real modes of the \ field separately, the full 
dissipation coefficient is well described by adding both contributions, in Eqs. (|4.9[) and (|4.15[) . as illustrated in Fig.fTUl 
for m x /T = 10. 

In the limit m a <C T, we may then write the following expression for the dissipation coefficient: 



T ^ 4 /m x 



4 



0.015/i 2 



( e m x /T _ 1)2 1 + Q 53 ^ 2 + q (I^)" 4 + Q.002^6 (!^)~ 6 \m x 



T 



(4.16) 



In Fig. [TT] we plot this expression by varying m x /T and keeping h fixed, which is more useful for practical im- 
plementations of warm inflation given that the former depends on the ratio <p/T, which is (slowly) varying during 
inflation. 

The curves shown in Fig. [TT] clearly illustrate the different regimes analyzed earlier, with the low- momentum (pole) 
contribution dominating at large (small) masses. It is also clear that for large values of h the contribution from virtual 
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m x /T=W, m o./T=0.01 




0.1 1 10 100 



FIG. 10: Comparison between the numerical results and the best fit expression obtained by adding the low-momentum and 
pole contributions to the loop integral for m x /T — 10 and m a /T = 0.01. 
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FIG. 11: Full dissipation coefficient as a function of m x /T for different values of the effective coupling h = h\/Ny, using the 
expression inferred from the numerical results. 

modes is different depending on whether their decay width is small or large compared to the mass, with the narrow 
width regime corresponding to the monotonic (T/m x ) 2 -bchavior at large masses. 

Given that during slow-roll inflation m x /T may not vary significantly, it is possible for the narrow-width form of 
the low-momentum dissipation coefficient in Eq. f|4. 1 1[) to hold during the full inflationary trajectory. In general, 
however, one expects transitions between the different regimes to take place over the 40 — 60 e-folds of inflationary 
expansion, and the dissipation coefficient may not have a monotonic behavior as illustrated in Fig. 111! Nevertheless, 
it should be emphasized that the contributions from real and virtual modes may actually have comparable values, 
although for different values of m x /T. 

V. HIGHER-ORDER CORRECTIONS 

In implementations of warm inflation based on the supcrpotential in Eq. (|2.1[) . one typically finds that large values 
of the constant > 10 s — 10 6 are required for a sufficiently long period of inflation (see e.g. [HI)- Although some 
of the assumptions used in these models may be relaxed, such as for example assuming that all MSSM degrees of 
freedom are present in the thermal bath during inflation, our previous discussion shows that large values of the 
effective coupling h = h\ZWy are generically required, unless an unnaturally large number of X-multiplets is coupled 
to the inflaton field. Note that, in this latter case, radiative corrections to the inflaton potential may still be kept 
under control by tuning the coupling g, which docs not modify the dissipation coefficient in the low- momentum and 
narrow width regime [see e.g. Eqs. (|4.7p and (|4.11|) ]. 
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It is then natural to investigate how higher-order corrections modify our earlier discussion and, in particular, the 
computation of the dissipation coefficient. Firstly, as we have discussed earlier, the condition m a < T should hold 
for large couplings and field multiplicities, and in fact the next-to-leading-order corrections come from the 'sunset' 
diagram shown in Fig. [5] and are proportional to h A Ny = h 2 h 2 , which can be made small in the limit of small couplings 
and large multiplicities such that h remains large. Similarly, the next-to-leading order corrections to the mass of the 
heavy fields, m x , arising from the bi-quadratic and quartic interactions of the <Ji fields, are Q(h 4 Ny , h 6 N Y ) and can 
also be made small in this limit. Moreover, for large tree-level masses, m x > T, thermal corrections are negligible up 
to large values of the effective coupling. 

The most relevant effect of higher-order corrections are, in fact, modifications to the decay width of the x-helds due 
to the self-interactions in the light scalar sector. These allow for emission of additional scalar pairs in the final decay 

state, with \ ~~ ^ <?iO'i(o'jO'j){Gk< J k) ■ • •> k, . . . = 1, • • • Ny, which will necessarily increase the total decay width. The 
effect of these higher-order many-body final states is hard to estimate in general due to the complexity of the resulting 
phase space, but to estimate this we have computed the leading corrections from four-body decays for on-shcll bosons 
at zero temperature. We describe this computation in detail in Appendix A, yielding a 'radiator' factor relative to 
the two-body partial decay width T x = RT^ of the form: 

r = i+ 1(1) 4n - + — m 

for Ny ^> 1. The leading correction is thus also proportional to h 4 Ny, so that it may be neglected for h <C 1 even 
though h > 1. 

On the other hand, for (h/ Air) 4 Ny > 1, we expect the radiator factor to become large, with significant contributions 
from an increasing number of many-body final states. Finite temperature effects may also increase this factor, as 
suggested by the leading order result in Eq. (|4.4[) . By enhancing the decay width, these corrections will also increase 
the dissipation coefficient in the narrow-width regime, as suggested by the results derived in the previous section. 
Although a full determination of the radiator factor is outside the scope of this work, we expect our results in this 
regime to be modified by replacing h 2 — > h 2 R, with the narrow width expressions now holding for m x /T > h\fR 1 
in which case the higher-order corrections to fh x have a negligible effect according to the above discussion. In fact, 
one should note that (h/Air) 4 Ny must not be too small or otherwise the thermalization rate of the light scalars will 
be suppressed, according to Eq. (|2.7[) . which suggests further exploitation of higher-order radiative corrections in the 
light scalar sector. 



VI. CONCLUSION 



In this work we have extended earlier analyses of dissipativc effects in supersymmetric warm inflation in the low- 
temperature regime, m x < T, to generic couplings and field multiplicities. As found in earlier analyses, the dissipation 
coefficient arising from the two-stage scalar interactions ^> — > % — > at differs significantly depending on whether the 
heavy fields \ are produced on- or off-shell. We have computed the dissipation coefficient numerically and obtained 
expressions that accurately describe our numerical results in both cases. Our full result, corresponding to a sum of 
both real and virtual contributions, is given by Eq. (|4.16[) in the limit m a <C T . We also find that the dissipation 
coefficient exhibits only a mild logarithmic dependence on the mass of the light fields, with m a < T even for large 
values of the coupling constant. 

We have found that, for a given value of the effective coupling h = h\/Ny, the production of on-shell x-modes from 
the rolling inflaton field is the dominant source of dissipation in the small mass limit, m x /T < 10, consistently with 
earlier analysis. For larger masses, the dominant process is the production of virtual x-pairs which then decay into 
the light degrees of freedom. In this case, the dissipation coefficient takes different forms depending on whether the 
field has a narrow or broad decay width, yielding T ~ C^T 3 /^ 2 and T ~ D^^/T 3 , respectively. Although in the 
narrow width regime we obtain Cf, cx h 4 for h < 1 and finite m a /T, in agreement with [HI, [35j], we find C$ oc h 2 
for large effective couplings, which is the regime where dissipative effects may be more pronounced. This constant is 
also enhanced by the multiplicity of the heavy fields, Nx, and by radiative corrections leading to many-body decay 
channels, yielding a radiator factor R that we have estimated to leading order. Taking all these considerations into 
account, we have: 

C ^ 4.7 RN x Ny , (6.1) 



with the narrow width approximation holding for m x /T > hy/NyR. Notice that in this regime the dissipation 
coefficient is independent of the coupling g between the inflaton and the heavy fields, which allows for large dissipative 
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effects, while keeping radiative and thermal corrections to the inflaton potential under control. Moreover, this shows 
that strong dissipative effects may be achieved by moderately large field multiplicities and couplings. In particular 
we expect radiative corrections to become significant for (/i/47r) 4 i\V > 1, which may in fact allow for warm inflation 
models with a small number of fields. 

Finally, one should emphasize that, despite the supersymmetric nature of the underlying model, interactions in 
the fermionic sector are inherently different from the scalar particles, as a consequence of supersymmetry breaking 
at finite temperature and energy density during warm inflation. In particular, both the dissipative production and 
the thcrmalization rate for light fermionic degrees of freedom are suppressed for m x /T > 1, so that fcrmions are a 
negligible (and most probably non-thermal) component of the radiation bath during inflation. We nevertheless expect 
fermions to be produced once the mass of x an( i V'x nas decreased sufficiently towards the end of inflation, which 
may happen in models where ip decreases or when there are additional tachyonic contributions to m x , as in hybrid 
inflation models. 

In this work, we have not considered the possible embedding of the supersymmetric model in Eq. <|2 . 1 1) in a gauge 
theory, mainly for simplicity, since the inclusion of gauge fields make the analysis much more involved. However, for 
minimally coupled x fields, we do not expect direct decays into gauge fields to be significant, since the only two-body 
process x ~~ ^ X7i although possible at finite temperature, is suppressed by the large x mass. Note that gauge bosons 
may nevertheless be emitted from the a scalars, if they are charged, or produced thermally, but this will necessarily 
involve a separate gauge coupling that can be made small enough for us to neglect these processes. 

We have shown in this work that by further exploiting regimes of parameters involving not too large coupling 
constants, such that perturbation theory remains valid, or at least that higher-order corrections are not large enough 
to spoil the lowest-order predictions, it is possible to enlarge the parameter space useful for warm inflation model 
building, by allowing for stronger dissipation at the cost of not too large field degeneracies, while thermal corrections 
to the inflaton potential arc still kept under control. Thus, we hope that this work will motivate further investigation 
of warm inflation model building and of the importance of dissipative effects in modern cosmological problems, both 
during and after the early inflationary phase. 
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Appendix A: Four-body decay at zero temperature 

The leading correction to the decay width of the heavy \ fields arises from the 4-body decays x < J i G \< J j' J ji 
illustrated by the diagram below. 



X.- -*- 



9 x 

-> c 

s 

s 

X 



FIG. 12: Feynman diagram for the 4-body decay x ~~ > <J i a \ (J j (J j- 

Here we will compute the associated partial decay width for an 'on-shell' heavy particle at T = 0. Assuming hi = h 
for simplicity, the corresponding amplitude is given by: 

M (4 > = ~(m x h 2 )h 2 „ 1 , , (Al) 
q z — m£ 
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where q 1 = (p — pi) is the momentum transfer. For the kinematics of the decay we follow [4l|, so that: 

( 4 ) X/3 12 /3 34 2|M (4 )| 2 
dr vto) = N fi s — ET^— ds 12 ds 34 d cos OkI cos 6 3 d(f> , A2 

where /?i2(34) is the velocity of particle 1 (3) in the CM frame of particles 1 and 2 (3 and 4), s\ 2 = (pi +P2) 2 , 
S34 = (p3 +P4) 2 , 6' 1 ( 3 ) is the angle between Pi(3) and Pi(3) + P2(4) and 4> is the angle formed by the planes (pi,P2) 
and (p3,P4). We have also defined: 



X = ^\jm\ + s\ 2 + s 2 4 - 2m 2 x (s 12 + s 34 ) - 2si 2 s 34 . (A3) 

The symmetry factor Sij accounts for the number of identical bosons in the final state, with 5 = 3! for i = j and 
5 = 2! for i =/= j, and we multiply the squared amplitude by a factor 2 to account for the fact that the (JjCJj pair may 
emitted from either primary cr_ field. 

For simplicity, we take the limit m a — > 0, in which we can show the following kinematical relations: 



Ei + E 2 = 



m 2 , + S12 - S34 
2m x 



sia = 2£ 1 (£i+.E 2 -|Pi2|cos0 1 ) 



IP12I = |pi+p 2 | = v / (-Ei+£ ; 2) 2 -si 2 . (A4) 



These then allow one to write: 



2 2 o P 2 ( m - sia - s 3 4 - 2m x cos gjjgia | ^ 

g = m - 2m x Ei = m — f- — . A5 

x x ym 2 , + S12 - S34 - 2m x cos6li|pi2|y 

Hence, = 9 2 (si2, S34, cos^i), so that the integrals in #3 and tfi are trivial. The remaining integrals can be performed 
by taking into account that < S\ 2 < to 2 , < S34 < (to x — ^/s 12 ) 2 and — 1 < cos^i < 1. 

The integrals actually exhibit logarithmic divergences in the limit m a — > 0, corresponding to the case where the 
secondary emitted scalars are soft and collinear with one of the primary scalars. As for example in QED and QCD, 
these soft collinear divergences will cancel against contributions from virtual scalars, i.e. when the OjOj pair is emitted 
and re-absorbed (see e.g. [HI), a result which should hold also at finite temperature [43[. Keeping only the finite 
terms and summing over all possible decay channels, we obtain: 

r(^A(3A y -2)(Ay r ( 2 ), (A6) 

(2) _____ 

where T x is given by the zero-temperature limit of Eq. (I4.3[) for an on-shcll \ boson. 
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